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Abstract 

Background: Over the last 30 years, extensive studies have revealed the crucial roles played by microbes in 
aquatic ecosystems. It has been shown that bacteria, viruses and protozoan grazers are dominant in terms of 
abundance and biomass. The frequent interactions between these microbiological compartments are responsible 
for strong trophic links from dissolved organic matter to higher trophic levels, via heterotrophic bacteria, which 
form the basis for the important biogeochemical roles of microbial food webs in aquatic ecosystems. To gain a 
better understanding of the interactions between bacteria, viruses and flagellates in lacustrine ecosystems, we 
investigated the effect of protistan bacterivory on bacterial abundance, production and structure [determined by 
16S rRNA PCR-DGGE], and viral abundance and activity of two lakes of contrasting trophic status. Four experiments 
were conducted in the oligotrophic Lake Annecy and the mesotrophic Lake Bourget over two seasons (early spring 
vs. summer) using a fractionation approach. In situ dark vs. light incubations were performed to consider the effects 
of the different treatments in the presence and absence of phototrophic activity. 

Results: The presence of grazers (i.e. < 5-Lim small eukaryotes) affected viral production positively in all 
experiments, and the stimulation of viral production (compared to the treatment with no eukaryotic predators) was 
more variable between lakes than between seasons, with the highest value having been recorded in the 
mesotrophic lake (+30%). Viral lysis and grazing activities acted additively to sustain high bacterial production in all 
experiments. Nevertheless, the stimulation of bacterial production was more variable between seasons than 
between lakes, with the highest values obtained in summer (+33.5% and +37.5% in Lakes Bourget and Annecy, 
respectively). The presence of both predators (nanoflagellates and viruses) did not seem to have a clear influence 
upon bacterial community structure according to the four experiments. 

Conclusions: Our results highlight the importance of a synergistic effect, i.e. the positive influence of grazers on 
viral activities in sustaining (directly and indirectly) bacterial production and affecting composition, in both 
oligotrophic and mesotrophic lakes. 

Keywords: Lakes microcosm, spring-summer variations, bacterial production, viral production, bacterial community 
structure, grazers 
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Background 

The heterotrophic bacterial community is the most 
important biological compartment involved in the trans- 
formation and mineralization of the organic matter in 
aquatic systems. It also constitutes a key source of prey 
for higher trophic levels, i.e. primarily flagellates, but 
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also ciliates and the metazooplankton [1,2]. Our concep- 
tual understanding of the role of heterotrophic bacteria 
in pelagic systems and in global biochemical cycles is 
closely linked to our understanding of how their growth 
rate, abundance, distribution and diversity are controlled 
[3-5]. 

Different biotic and abiotic factors have been identi- 
fied as players acting on the activity and composition of 
the bacterial community, and resources (organic matter 
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and nutrients) are considered one of the main factors 
controlling this community [2,6]. However, the roles of 
bacterivory and viral lysis are not insignificant, and may 
also strongly affect bacterial abundance, activity and 
structure. Both heterotrophic nanoflagellate (HNF) graz- 
ing and viral lysis are known to be variable causes of 
bacterial mortality, and can be responsible for 10 to 60% 
of daily bacterial loss in lacustrine systems [e.g. [7]]. In 
addition, both processes can impact the size distribution 
of bacterial communities through 'size-selective mortal- 
ity' for flagellates [8,9] and 'host-specificity' for viruses 
[10]. Moreover, viruses can act indirectly on bacterial 
structure throughout the release of cell debris during 
lysis activity (enriching the pool of dissolved and parti- 
culate organic matter (DOM and POM) and inorganic 
nutrients) enhancing in fine growth and production of 
some bacterial groups [11,12]. Indeed, whether cells are 
grazed or lysed can have different ecological and biogeo- 
chemical consequences, as the implications for the mat- 
ter and energy flow through the microbial web will be 
very different [13,14]. Typically, high rates of viral cell 
lysis may generate a recycling of nutrients and organic 
matter at the base of the food web and therefore, less 
carbon and nutrients may reach higher trophic levels, a 
process referred to as the viral shunt [13,14]. In con- 
trast, if bacteria are grazed by flagellates, nutrients and 
energy can reach higher trophic levels via the connec- 
tion between the microbial loop and the classical food 
chain [15]. Thus, these processes can significantly influ- 
ence the production of dissolved organic carbon and the 
recycling of nutrients [14,16] and can impact/modify 
not only bacterial diversity [9,17] but also the relation- 
ship between diversity and ecosystem functioning [18]. 

A few studies have investigated the individual effects 
of flagellates or viruses on bacterial communities in 
terms of abundance, production and diversity (e.g. 
[7,10,19,20]). However, their combined effects on bac- 
teria, and the comparison between individual and com- 
bined effects are still limited [18,21,22]. According to 
these studies, both viral lysis and protistan bacterivory 
may act additively to reduce bacterial production and 
sustain diversity, which could explain the less pro- 
nounced blooming species in heterotrophic bacterio- 
plankton than in phytoplankton [22]. However, the 
opposite effect has also been reported [23]. Moreover, 
comparisons of the combined effects of viruses and fla- 
gellates on the bacterial community according to the 
trophic status of aquatic systems are scarce and until 
now, no information has been made available for lacus- 
trine systems. To the best of our knowledge, Zhang et 
al. [22] are the only authors who have investigated these 
effects taking into account a trophic range within a 
coastal ecosystem, and the same trend was highlighted 
[22]. According to these authors, a shift of predator 



control mechanisms from flagellates in oligotrophic sys- 
tems to viruses in eutrophic systems could explain the 
results. 

In this study, we collected samples from two peri- 
alpine lakes (Annecy and Bourget) with substantial dif- 
ferences in their trophic state (oligo- vs. mesotrophic, 
respectively) and we developed treatments with either 
individual or combined predators of the bacterial com- 
munity using a fractionation approach (i.e. a physical 
separation of virus-bacteria and the small eukaryotes). 
Our main goal was to examine the separated and com- 
bined effect of viruses, grazers and small autotrophs (< 
5 (im) on the bacterial abundance, production and 
structure, and to compare it in different environmental 
conditions. Since the importance of both predators (fla- 
gellates and viruses) as potential controlling forces of 
the bacterial community may display seasonal variations 
in these lakes [7,8,24], this study was carried out at two 
contrasting periods (early-spring vs. summer), character- 
ized by substantial differences in both the dynamics and 
structure of microbial communities and environmental 
conditions [8,25]. 

Our main findings are that both viral lysis and flagel- 
lated bacterivory act additively to sustain bacterial pro- 
duction, probably through a cascading effect from 
grazer-mediated resource enrichment, whereas their 
effects on the bacterial community structure remain 
more subtle. On the whole, the combined effects of 
viruses and flagellates showed the same trend in both 
lakes Annecy and Bourget. 

Results 

Initial conditions 

In situ characteristics of the study sites 

Lake Bourget is an elongated and north-south oriented 
lake situated in the western edge of the Alps (length 18 
km; width 3.5 km; area 44 km 2 ; volume 3.5 x 10 9 m 3 ; 
altitude 231 m; maximum depth 147 m; mean depth 80 
m; residence time 8.5 years). Lake Annecy is located in 
the eastern part of France, at a distance of approx. 50 
km from the former, (length 14.6 km; width 3.2 km; 
area 28 km 2 ; volume 1.2 x 10 9 m 3 ; altitude 447 m; maxi- 
mum depth of 65 m; mean depth 41 m; residence time 
3.8 years). From the end of March to mid-July (i.e. peri- 
ods during which experiments were conducted), in situ 
temperatures of the two study sites varied between 6.2° 
C and 20.4°C, while the dissolved oxygen varied more 
modestly, between 9.7 and 11.7 mg l" 1 (Table 1). Differ- 
ences in the concentration of nutrients (N0 3 , NH 4 and 
Ptot) between Lake Annecy and Lake Bourget were 
principally recorded during the early spring experiments 
(LAI and LB1, respectively), with values twice to three- 
times higher in Lake Bourget (LB1) than in Lake Annecy 
(LAI) (Table 1). Chi a concentration was relatively low 
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Table 1 Physicochemical and biological characteristics of 
the sampling sites (2 m depth) 



Parameters 




LAI 


LA2 


LB1 


LB2 


Sampling 




26/03/ 


10/07/ 


02/04/ 


1 7/07/ 


date 




2007 


2007 


2007 


2007 


Temperature 


°r 
L 


6.2 


1 9.6 


7.5 


20.4 


DO 


mg I" 1 


10.5 


9.7 


11.7 


10 


TOC 


mg I" 1 


1.7 


2.2 


2.1 


2.5 


N0 3 


mg I" 1 


0.2 


0.1 


0.5 


0.2 


NH 4 


M9 I" 1 


2.0 


1.0 


6.0 


4.0 


P0 4 


M9 I" 1 


2.0 


3.0 


4.0 


2.0 


r loiai 


i in r 1 


/.u 


o.u 


1 1 n 
z I .U 


o.u 


Chk? 


M9 I" 1 


0.7 


2.7 


1.2 


0.7 


Cyanobacteria 


10 4 cell 


9.0 ± 0.5 


15.0 ± 1.1 


2.0 ± 0.1 


12.0 ± 0.8 




ml" 1 










Het. Bacteria 


10 5 cell 


24.4 ± 0.3 


12.3 ± 0.4 


35.0 ± 1.2 


25.2 ± 2.0 




ml" 1 










Viruses 


10 7 part 


3.7 ± 0.1 


5.1 ± 0.4 


8.3 ± 0.3 


15.3 ± 0.7 




ml" 1 










HNF 


10 2 cell 


7.5 ± 1 .3 


6.9 ± 0.6 


2.6 ± 1 .3 


3.9 ± 1.5 




ml" 1 










PNF 


10 2 cell 


4.9 ± 1 .3 


18.0 ± 3.1 


1 .4 ± 0.2 


2.9 ± 0.5 




ml" 1 











DO, dissolved oxygen; Chi a, Chlorophyll a; TOC, total organic carbon; N0 3 , 
nitrate; NH 4 , ammonium; P total, total phosphorus; Het. Bacteria, heterotrophic 
bacteria; HNF, heterotrophic nanoflagellates; PNF, pigmented nanoflagellates. 
LA1, March sampling in Lake Annecy; LA2, July sampling in Lake Annecy; LB1, 
April sampling in Lake Bourget; LB2, July sampling in Lake Bourget. Values are 
means ± standard deviation of results from triplicate measurements. 

(i.e. < 2.8 (ig l" 1 ) for the four experiments (LAI, LA2, 
LB1 and LB2). The abundance of heterotrophic bacteria 
varied between 1.2 and 3.5 x 10 6 cell ml" 1 , viruses 
between 3.7 and 15 x 10 7 virus ml" 1 , heterotrophic 
nanoflagellates (HNF) between 2.6 and 7.6 x 10 2 cell 
ml" 1 , pigmented nanoflagellates (PNF) between 1.4 and 
18 x 10 2 cell ml" 1 , and picocyanobacteria between 2 and 
15 x 10 4 cell ml" 1 . These parameters were significantly 
different (ANOVA, P < 0.05, n = 12) between the four 
experiments (LAI, LA2, LB1 and LB2), indicating dis- 
tinct biological characteristics at initial sampling. Seaso- 
nal difference in the picocyanobacterial abundance was 
monitored (ANOVA, P < 0.05, n = 6) in both lakes 
(Annecy vs. Bourget), with values 1.6- to two-times 
higher in summer (LA2 and LB2) than in early spring 
(LAI and LB1). HNF and PNF abundances were signifi- 
cantly higher in Lake Annecy than in Lake Bourget 
(ANOVA, P < 0.05, n = 12, Table 1). In contrast, viral 
abundances were always lower in the oligotrophic Lake 
Annecy. 

Conditions in experimental bottles - Effect of filtration 

The < 5-|im prefiltration removed a relatively small frac- 
tion of both HNF and PNF (less than 20%), whereas the 
< 1.6-(im filtration removed, as expected, all of them 



(Table 2). At the start of the experiments, in VF 
(Viruses+Bacteria+Flagellates) and VFA (Viruses+Bac- 
teria+Flagellates+Autotrophs) treatments, HNF abun- 
dances varied between 2.5 x 10 2 cell ml" 1 (LB) and 6.5 x 
10 2 cell ml" 1 (LA), PNF between 1.1 x 10 2 cell ml" 1 (LB) 
and 14.4 x 10 2 cell ml" 1 (LA), and picocyanobacteria 
between 0.7 x 10 4 cell ml" 1 (LB) and 11.2 x 10 4 cell ml" 1 
(LA) corresponding to 52 to 72% of in situ abundances. 
Comparatively, a small fraction of the picocyanobacterial 
community passed through the < 1.6- (im filter and only 
0.1 and 0.8 x 10 4 cell ml" 1 were recorded in treatment 
V (only bacteria and viruses), i.e. 1 to 5% of in situ 
abundance (Tables 1 and 2). In contrast, filtration 
through 1.6 \im resulted in a small loss of bacterial and 
viral abundances (less than 14% and 20%, respectively) 
whereas after 5-(im filtration, loss never exceeded 4% 
for heterotrophic bacteria and 13% for viruses. At the 
beginning of the incubation, heterotrophic bacteria and 
viral abundances, in the four treatments of all experi- 
ments varied between 9.4 x 10 5 and 33.5 x 10 5 cell ml" 1 
and between 2.9 x 10 7 and 13.4 x 10 7 virus ml" 1 , respec- 
tively (Figure 1). Overall, we succeeded in obtaining 
incubations with strongly contrasting predator pressure 
(with or without) and, with negligible loss to the abun- 
dances of both bacteria and viruses, when compared to 
in situ conditions. 

Effect of treatments on viral abundance and production 

Viral abundance only varied by a small degree (between 
2.9 x 10 7 and 4.6 x 10 7 virus ml" 1 ) in Lake Annecy, 
while it varied greatly in Lake Bourget particularly dur- 
ing the LB2 experiment (Figure 1). In both LAI and 
LA2 experiments, the temporal trend of viral abundance 
revealed different patterns according to the treatment: 
viral abundance increased in VF and V treatment, while 
in the VFA treatment no significant evolution (ANOVA, 
P > 0.05, n = 9) was recorded (Figure 1). In Lake Bour- 
get, viral abundance increased during the four days of 
incubation in all treatments, except in treatment V of 
the LB1 experiment. At the end of incubation, the 
increase in viral abundance in VF and VFA was signifi- 
cantly higher than in treatment V (ANOVA, P < 0.01, n 
= 9) in LAI (+39% and +16%, respectively), LB1 (+34% 
and +27%, respectively) and LB2 (+47% and +61%, 
respectively) (Figure 2D). However, the opposite was 
true for LA2 (-6%, ANOVA, P < 0.05, n = 9). The sti- 
mulation of viral abundance was 3-fold higher in Lake 
Bourget (average +29%) than in Lake Annecy (average 
+8%) (t test, n = 24, P < 0.001) and significantly differ- 
ent between seasons for each lake (t test, n = 12, 
p < 0.01). 

Similarly to viral abundance, viral production 
increased without exception from the beginning to the 
end of the experiments, particularly in VFA and VF 
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Table 2 Picocyanobacteria (Synechococcus spp), pigmented and heterotrophic nanoflagellates (PNF and HNF) 
abundances at the beginning (t 0 ) and at the end (t fina i) of each incubation in the different experiments 



Experiments and treatments 






Abundance of: 








Picocyanobacteria (10 


4 cell ml 1 ) 


PNF (10 2 cell ml' 1 ) 


HNF (10 2 


cell ml' 1 ) 


to 


tfinai 




Minal 


to 


tfinai 


LA1 


















V 


0.1 ± 0.0 


0.3 ± 0.0 


U.U ± U.U 


U.U ± U.U 


0.0 ± 0.0 


0.0 ± 0.0 




VFA 


6.5 ± 0.1 


7.5 ± 0.1 


/I C -u 1 Q 


4.o ± U.j 


6.2 ± 1.3 


8.1 ± 1.4 




VF 


5.5 ± 0.1 


2.4 ± 0.2 


4.Z x U.z 


c. c. _i_ n a 

0.0 x U/+ 


6.5 ± 0.9 


8.0 ± 2.6 


LA2 


















V 


0.8 ± 0.4 


0.3 ± 0.2 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 




VFA 


10.2 ± 0.1 


15.8 ± 0.1 


14.4 ± 0.6 


28.5 ± 1 .3 


5.6 ± 0.2 


11.1 ± 0.8 




VF 


1 1 .2 ± 0.4 


6.3 ± 0.3 


14.0 ± 0.4 


19.1 ± 0.1 


5.4 ± 0.3 


13.5 ± 0.8 


LB1 


















V 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 




VFA 


0.8 ± 0.0 


1.5 ± 0.1 


1 .3 ± 0.5 


8.7 ± 0.5 


2.5 ± 0.5 


12.0 ± 1.7 




VF 


0.7 ± 0.2 


0.4 ± 0.3 


1.1 ± 0.7 


6.5 ± 0.2 


2.9 ± 0.6 


12.4 ± 0.2 


LB2 


















V 


0.3 ± 0.0 


0.5 ± 0.1 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 


0.0 ± 0.0 




VFA 


7.3 ± 0.1 


16.6 ± 2.1 


2.5 ± 2.8 


7.5 ± 8.9 


3.6 ± 4.1 


20.7 ±11.7 




VF 


7.1 ± 0.7 


3.1 ± 0.2 


3.1 ± 1.5 


12.5 ± 0.9 


3.9 ± 4.0 


13.8 ± 9.0 



V, Viruses+Bacteria treatments; VFA, Viruses+Bacteria+Flagellates+Autotrophs treatments; VF, Viruses+Bacteria+Flagellates treatments. LA1, LA2, LB1, LB2: 
abbreviations as in Table 1. 



treatments whatever the lake (ANOVA, P < 0.05, n = 9). 
Viral production varied between a minimum of 3.2 x 
10 5 virus ml" 1 h" 1 (LA2) and a maximum of 4.7 x 10 6 
virus ml" 1 h" 1 (LB2) (Figure 3), which corresponded to 
3.5 x 10 5 and 47.4 x 10 5 cells lysed ml" 1 d" 1 , respectively 
(Table 3). Viral production in VFA and VF were, in 
most cases, significantly different (ANOVA, P < 0.001, n 
= 18), over the course of the incubation, being on aver- 
age 21% higher (range 7-53%) than in V treatments in 
both lakes (Figure 2E). Stimulation of viral production 
seemed to be significantly higher (t test, P < 0.0001, n = 
24) in Lake Bourget (average +30%) than in Lake 
Annecy (average +11%), while no significant seasonal 
differences (t test, P > 0.05, n = 12) were recorded for 
either lake. 

Effects of treatments on bacterial abundance, production 
and mortality 

Bacterial abundance increased throughout the experi- 
ments, particularly during the LB2 experiment (Figure 
1). Concentrations were significantly higher in VFA and 
VF than in treatment V (ANOVA, P < 0.05, n = 18). 
Concentrations in VFA and VF were in most cases simi- 
lar in Lake Annecy, when compared to each other 
(ANOVA, P > 0.05, n = 18), in contrast to the signifi- 
cant differences observed in the samples issued from 
Lake Bourget, with higher bacterial abundance in treat- 
ment VFA than VF. At the end of the incubation, the 
increase in bacterial abundance (comparison of 



treatments V and both VF and VFA between day 0 and 
day 4) in treatment VFA was significantly higher than in 
treatment V (ANOVA, P < 0.01, n = 9) (Figure 2A). In 
the four experiments, bacterial abundance was signifi- 
cantly higher (by up to 9% to 53%) (t test, P < 0.05) in 
treatment VFA than in V. In the VF treatment, bacterial 
abundance was significantly higher (t test, P < 0.05) in 
LA2 (up to 35%), LB1 (up to 30%) and LB2 (up to 19%) 
than in treatment V. No significant difference was 
observed in LAI (t test, P>0.8). Stimulation of bacterial 
abundance was significantly different between lakes (t 
test, P < 0.001, n = 24) (+38% in Lake Bourget and 
+ 14% in Lake Annecy) and between seasons with high- 
est values measured in summer (+59% in Lake Bourget 
and +26% in Lake Annecy). 

During the incubation period, bacterial production 
fluctuated between 0.5 and 0.9 ugC l" 1 h" 1 in LAI, 0.8 
and 2.3 ugC l" 1 h" 1 in LA2, 1.2 and 3.1 ugC l" 1 h" 1 in 
LB1 and between 3.2 and 7.8 ugC l" 1 h" 1 in LB2 (Figure 
3). Following bacterial abundance evolution, a significant 
increase in the bacterial production (ANOVA, P > 0.05, 
n = 27) was also recorded throughout the period of 
incubation. For both lakes, bacterial production was 
often higher in treatment V than in both VFA and VF 
during the early spring experiments (LAI and LB1). 
After 96 h of incubation, the stimulation of bacterial 
production (comparison of variation of the viruses treat- 
ment (V) and the grazers treatments (VFA and VF)) was 
observed in all experiments and averaged 27% in 
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0 48 

Time (hours) 




0 48 

Time (hours) 



v □ WA □ w 

Figure 1 Time-course of viral abundance (10 7 virus ml 1 ) and bacterial abundance (10 6 cell ml 1 ) in the four experiments during the 
incubation period. Values are given as mean ± standard deviation of triplicate incubations. Asterisks indicate sampling time point for which the 
VFA and VF treatments were not significantly different from the V treatment (ANOVA, P > 0.05, n = 9). Note that the panels have different scales. 
LA1, LA2, LB1, LB2: abbreviations as in Table 1. 
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3 



0 

E 

0 



o 

3 



CO 



10a 
so 

60- 
40" 
20- 
0 

-20- 
-40- 
-60 



A 



"Dp" 



f 



On 



ffi. 



EL 





LA1 LA2 LB1 LB2 
48h 



jo 



LA1 ' l3\2 ' LB1 LB2 
96h 



□ Treatment Viruses+Flagellates+Autotroph (VFA) 

□ Treatment Viruses+Flagellates (VF) 



100 

30- 
60" 
40" 
20 
0 

-20 

-4a 

-60 

□ Treatment Viruses+Flagellates+Autotroph (VFA) 

□ Treatment Viruses+Flagellates (VF) 

Figure 2 Stimulation of bacterial abundance (A) 7 production (B) and lysis mortality (C) 7 and viral abundance (D) and production (E) 

Values are given as mean ± standard deviation of triplicats incubations from VF and VFA treatments and average values from the V treatment. 
ANOVA were preformed between treatment V and the two other treatments (VF and VFA), and the comparison without significance (P > 0.05) 
was indicated with asterisks. 



treatment VFA and 20.8% in treatment VF when com- 
pared to V (Figure 2B). The highest stimulation was 
observed in VFA during the LB2 experiment (51%). 
Overall, the bacterial production was significantly differ- 
ent (ANOVA, P < 0.001, n = 27) between the three 
treatments for the four experiments, with the highest 
values observed in most cases in VFA and VF (Figures 2 
and 3). In contrast to the bacterial abundance, a signifi- 
cant difference in the stimulation of bacterial production 
was only noted between seasons (t test, P < 0.001, n = 
12), with the highest values for summer experiments 



(+33.5% and +37.5% for Lake Bourget and Lake Annecy, 
respectively). 

Bacterial growth rate fluctuated between 0.1 and 0.7 
d" 1 after either 48 h or 96 h of incubation (Table 3), 
with the lowest values recorded during early spring 
experiments (LAI and LB1). The presence of flagellates 
did not induce a reduction of bacterial abundance and 
the estimation of bacterial loss rates over time gener- 
ally led to negative values, showing enhanced bacterial 
growth. In Lake Annecy, this positive impact on bac- 
terial growth was only significant in the LA2 
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LA1 



* * 



l 



I 




E 

3 



O 



O 8 
-♦— ' 
O 
3 

"8 - 



> 



2 



LA2 



LB1 



f 



IT 1 




o 1 



c 

.2 o 



LA2 




60 



LB2 



i 



ft 



0 48 

Time (hours) 



96 




VZA V □ VFA E3 VF 



0 48 

Time (hours) 



Figure 3 Time-course of viral production (10 virus ml h ) and bacterial production (ugC I h ) in the four experiments during the 
incubation period. Values are given as mean ± standard deviation of triplicate incubations. Asterisks indicate sampling time point for which 
VFA and VF treatments were not significantly different from the V treatment (P > 0.05, n = 9, ANOVA). Note that the panels have different scales. 
LA1, LA2, LB1, LB2: abbreviations as in Table 1. 



Berdjeb et al. BMC Microbiology 201 1, 1 1:88 
http://www.biomedcentral.eom/1 471 -21 80/1 1 /88 



Page 8 of 18 



Table 3 Bacterial growth rate (r), loss rate, virus-induced mortality and lysis activity rates after 48 h and 96 h 



Experiment/Treatment Growth rate (r) (d 1 ) Loss rate of bacteria (d 1 ) Lysis mortality (d 1 ) Lysis rate activity (10 5 cell ml" 1 d" 1 ) 





48 h 


96 


h 


48 h 


96 h 


48 


h 


96 h 


48 


h 


96 


h 


LA1 




































V rA 


0.12 


± 0.05 


n 1 a + 

U. I JZ 


U.U I 


-0.03 


± 0.09 


-0.06 


± 0.04 


0 1 9 4- 
U. I o X 


o m 

U.U I 


0.21 


± 0.01 


D.yU X 


0 1 A 
U. I o 


4.60 


± 


0 OA 
U.U^f 


\/F 
v r 


0.09 


± 0.06 


n 1 n + 


u.uo 


0.01 
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The significant difference between bacterial growth rate in V treatment and VFA/VF treatments was tested using ANOVA. *, P < 0.05; **, P < 0.001. 



experiment (ANOVA, P < 0.05, n = 6), and was 
observed in both VF (-0.1 d" 1 ) and VFA (-0.1 d" 1 ). In 
Lake Bourget, the two experiments (LB1 and LB2) 
showed the same effect on bacterial growth, with the 
highest values observed in VFA treatment (-0.2 d" 1 , 
ANOVA, P < 0.001, n = 6). 

Bacterial mortality due to viral lysis activity was esti- 
mated to range between 0.2 d" 1 and 2.2 d" 1 (Table 3) 
with the highest values obtained during summer experi- 
ments (LA2 and LB2). Differences between V and VFA/ 
VF treatments indicated a significant increase in the 
lysis mortality rate after 48 h incubation in both LB1 
(+28%) and LB2 (+43%) and this enhancement was 
maintained until the end (96 h) (Figure 2C). In LAI and 
LA2, a significant difference between V and the other 
treatments was observed at the end of incubation, 
accompanied with an increase in lysis mortality rate in 
LAI (+11%), and a decrease in LA2 (-7%). 

Effects of treatments on the bacterial community 
structure 

Figure 4 shows the PCR-DGGE patterns of the bacterial 
community structure at the start and end of incubation 
for the three treatments and the four experiments. 
Between 17 and 26 bands were found in treatment V, 
between 18 and 28 in VF and between 18 and 27 in VFA 
(Figure 4 and Table 4). The number of common bands 
found in the three treatments for each experiment 



represented between 24 and 49% (average 40.5%, Table 4). 
Between 0 and 3 bands (average 3.8%) per experiment 
were specific to V. Between 0 and 2 bands (average 2.3%) 
and between 1 and 4 (average 6.5%) bands were specific to 
VF and VFA, respectively (Table 4). 

Among the eight sequenced bands, B2 (excised from 
the LAI experiment) was present in all treatments. 
From this band, ten sequences out of 12 obtained were 
related to the genus Curvibacter (class of P-proteobac- 
teria), the two other sequences corresponding to the 
genus Burkholderidia (class of P-proteobacteria) (Table 
5). Three other sequenced bands were visible in all 
treatments but they increased significantly in intensity 
at the end of incubation (both B3 and B4 in V fina i of 
LAI, B8 in VF fina i of LB2). These three excised bands 
were related to the phylum Actinobacteria (with B3 
affiliated to the clade acl) (Figure 4 and Table 5). 
Finally, the three last bands chosen to be sequenced 
appeared (B5 in V fina i and VF final of LA2) or disappeared 
(both B6 and B7 in VFA fina i of LB1) at the end of incu- 
bation (Figure 4). These ones were all affiliated to the 
phylum Actinobacteria (as were 85% of the sequenced 
DGGE bands). Note that the excised band Bl (LAI 
experiment), related to the phylum Cyanobacteria 
(Table 5), disappeared at the end of the incubation in 
both VF and V treatments. 

Cluster analyses based on quantification of the band 
position and intensity (Figure 5) showed that, for each 



Berdjeb et al. BMC Microbiology 201 1, 1 1:88 
http://www.biomedcentral.eom/1 471 -21 80/1 1 /88 



Page 9 of 18 



LA1 LA2 LB1 LB2 

VFA flnal VF fral V 0 V^f VF 0 VFA final VF fmai V 0 VFA, in „ VF„„„ V n VFA-„,VF-„, \ 




Figure 4 Bacterial community structure at the beginning 
(referred to as '0') and at the end (96 h, referred as 'final') of 
the incubation, visualized by DGGE of PCR-amplified partial 
16S rRNA genes, and the position of the different bands 
excised and sequenced. (B1 to B8, see Table 5). V 0 and V Fjna |, 
treatment Viruses+Bacteria at the beginning and the end of 
experiments; VF 0 and VF fina |, treatment Viruses+Bacteria+Flagellates 
at the beginning and the end of experiments, VFA fjna |, treatment 
Viruses+Bacteria+Flagellates+Autotrophs at the end of experiments 



lake, the bacterial community structure was clearly dif- 
ferent according to the period (early spring/summer) 
(Figure 5). In lake Annecy, different responses could be 
observed following our manipulations: In spring, the 
effect of time (separation between tO and tf, whatever 
the planktonic size fraction) clearly affected the bacterial 
structure. In summer, the differences in the structure 
induced by the size fractionation were the strongest, and 
sample discrimination was clearly linked to the fractio- 
nation (1.6 vs. 5 urn). Similar patterns were obtained for 
Lake Bourget in summer. Finally, treatment VFA was 
highly divergent from V and VF (between 42% and 58% 
of similarity) during the early spring experiment for 
Lake Bourget. 



Table 4 Number and specificity of bacterial DGGE bands 
in the four experiments 





LAI 


LA2 


LB1 


LB2 


Total Bands 


37 


31 


34 


31 


Treatment Viruses 


22-23 


17-20 


26 


17-25 


Treatment Viruses + Flagellates 


25-28 


18-19 


27 


19-23 


Treatment Viruses + Flagellates 
+Autotrophs 


25-28 


18 


24-27 


19-27 


Bands common in all treatments 


47% 
(16) 


42% 
(13) 


49% 
(18) 


24% 
(13) 


Bands specific to: 










Treatment Viruses 


0 


3% (1) 


9% (3) 


3% (1) 


Treatment Viruses + Flagellates 


0 


6% (2) 


3% (1) 


0 


Treatment Viruses + Flagellates + 
Autotrophs 


11% 
(4) 


6% (2) 


6% (2) 


3% (1) 


Two combined treatments (VFA + 
VF) 


8% (3) 


6% (2) 


15% 
(5) 


23% 
(7) 



Discussion 

Experimental approach 

In order to study the influence of both predation pres- 
sure and the autotrophic activity on bacterial commu- 
nity of Lakes Annecy and Bourget, we carried out a 
fractioning approach and performed incubation in either 
darkness or ambient light. The originality and strength 
of this study comes from the fact that such experiments 
were conducted (i) in two ecosystems with either oligo- 
trophic or mesotrophic status and (ii) at two distinct 
periods of the year (i.e. early-spring and summer) where 
microbial planktonic dynamics and composition are 
likely to display clear differences [8,24,25]. Although the 
use of microcosms may introduce some bias into the 
development of microbial communities compared with 
those occurring naturally in the field (due to 

Table 5 Phylogenetic information about the OTUs 
corresponding to the excised and sequenced DGGE 
bands 

Bands Number of OTUs Nearest uncultivated 

N° sequenced species accession no°,% 

clones similarity 



B1 



B2 



B2 



B2 



B3 



B4 



B5 



B5 



B7 



12 



10 



10 
1 

12 



Phylum: 
Picocyanobacteria 
Synechococcus sp 

Class: p- 
proteobacteria 
Genus: 
Curvibacter 

Class: p- 
proteobacteria 

Genus: 
Burkholderia 

Class: p- 
proteobacteria 

Genus: 
Burkholderia 

Phylum: 
Actinobacteria 
Clade: acl 

Phylum: 
Actinobacteria 
Unidentified 

Phylum: 
Actinobacteria 
Unidentified 

Unclassified 
bacteria 

Phylum: 
Actinobacteria 
Unidentified 

Phylum: 
Actinobacteria 
Unidentified 

Phylum: 
Actinobacteria 
Unidentified 

Unclassified 
bacteria 



AY224199, 98% 



EU703347, 98 
EU642369, 99% 



EU642141, 98% 



EU801155, 97% 
EU63973669, 96% 



FJ916243, 99% 



FN668296, 99% 



FN668268, 100% 



FJ916291, 99% 



DQ3 16369, 99% 



AJ575506, 99% 



Berdjeb et al. BMC Microbiology 201 1, 1 1:88 
http://www.biomedcentral.eom/1 471 -21 80/1 1 /88 



Page 10 of 18 



0.2 



0.4 



0.6 



0.8 



1 



B 



0.2 



0.4 



0.6 



0.8 



-V 0 
-VF 0 

■ VFAfjnai 

" Vfjnal 

- VF fina | 

-v 0 

- Vfinal 

-VF 0 

-VFA final 

-VF final 



LA1 



y LA2 



1 



• VFA fina | 

■V 0 

' Vfi na | 
VFfina, 

VF 0 

■V 0 
•VF 0 

■ Vfinal 

• VFA fin a| 

VF fina | 



V» LB1 



LB2 



Figure 5 Cluster analysis of DGGE profiles based on band position and intensity. Scale bars indicated the Bray-Curtis similarity index 
in Lake Annecy (A) and Lake Bourget (B). V 0 and V Fina |, treatment Viruses+Bacteria at the beginning and the end of experiments; VF 0 and 
VFf jna | ; treatment Viruses+Bacteria+Flagellates at the beginning and the end of experiments, VFA fjna i, treatment Viruses+Bacteria+Flagellates 
+Autotrophs at the end of experiments. 



confinement and handling effects), these experimental 
tools are still very useful for investigating how processes 
such as mortality factors induce temporal variation in 
bacterial dynamics, structure and activity [26]. Incuba- 
tion time (4 days) coupled with the volume of micro- 
cosms (2.5 L) considered in this study have previously 
been used successfully in other experimental studies 
[18,22]. We assumed that our design was thus realistic 
enough compared to the generation time of microorgan- 
isms and aimed to obtain significant changes in bacterial 
and viral activity [27] . A comparison of virus and flagel- 
late abundances at the onset of the experiments with in 
situ conditions and among treatments with different 
viral and flagellate effects was successful. However, the 



experimental protocol resulted in a reduction of HNF at 
the start of the experiment and we thus might have 
underestimated their influence. Clear effects of HNF 
were observed at the end of the experiment, when fla- 
gellate abundance was about twice as high as in situ 
(Tables 1 and 2). 

Grazing effect on viral activity 

According to the model of Miki and Yamamura [28], 
grazers should reduce the role of the viral loop. As viral 
production was higher in VF and VFA treatments than 
in V in most cases, particularly after 48 h incubation, 
our results instead suggested a stimulation of viral activ- 
ity in the presence of grazers. Such a stimulation of viral 
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production by the presence of small eukaryotes (grazers) 
was observed in all experiments for the two lakes. These 
results corroborate the findings of Jacquet et al. [27] 
who observed a clear and positive relationship between 
flagellate concentration and VIBM (virus-induced bac- 
terial mortality) in Lake Bourget (r = 0.99, P < 0.05) at 
three different periods of the year (winter, spring and 
summer), suggesting a synergistic cooperation between 
grazer and virus activity. Our new results extend the 
occurrence of this process at other periods of the year 
and in the oligotrophic Lake Annecy. Similar beneficial 
effects of protozoan grazing on viruses have been 
reported in various lacustrine systems with different 
trophic statuses [21,23,26]. This means that the trophic 
status cannot be 'used' as an environmental factor to 
change the balance between positive and negative effects 
of flagellates on viruses [29], and it is likely that there 
are probably different processes involved in enhancing 
viral activities in response to grazing activity [21]. 

To the best of our knowledge, Simek et al. [19] were 
first to suggest that protozoan grazing may influence 
and increase viral lysis. From that time, other studies 
reported such a synergistic effect in contrast to fresh- 
water systems [21,26,27]. Nevertheless, an antagonistic 
interaction between these two compartments was also 
noted elsewhere [30,31]. Mechanisms by which HNF 
affect viral activity are still unclear and many hypotheses 
have been proposed to explain such a cooperative inter- 
action (reviewed by Miki and Jacquet [29]). In brief, 
grazing activity could stimulate bacterial growth rates, 
by releasing organic and inorganic nutrients. Higher 
bacterial growth rates might be associated with 
enhanced receptor formation on cell surface which may 
result in a greater chance of phage attachment and in 
fine higher infection frequencies. Thus, grazer stimula- 
tion of viral proliferation could occur through cascading 
effects from grazer-mediated resource enrichment [23]. 
We observed, in this study, a strong stimulation of bac- 
terial production in treatments with grazers which may 
corroborate this assumption in both lakes. A link 
between infection and host production has been 
reported previously (summarized in Weinbauer [11]) 
and, recently, experimental studies showed that viruses 
may preferentially lyse active cells [18,32]. Our results 
showed that autotrophic activity contributed to this sti- 
mulation, mainly in the early summer experiment (for 
both lakes), while heterotrophic flagellates were always 
involved in this positive feedback. A shift in the bacterial 
community structure could also contribute to the syner- 
gistic interaction observed in this study. According to 
Weinbauer et al. [33], grazing can favour species which 
are capable of rapid growth and resistant to grazing, 
which may stimulate viral infection. Our results also 
seemed to support this hypothesis since both high 



bacterial production and specific bands were only 
observed in treatments VF and VFA. 

Stimulation of viral production was much more vari- 
able between lakes than between seasons and it was 
clearly higher in Lake Bourget. This suggests that envir- 
onmental conditions encountered in the mesotrophic 
system might promote higher viral activity compared to 
more oligotrophic conditions. This hypothesis agrees 
with Lymer et al. [34] or Pradeep and Sime-Ngando 
[26] who observed, during a microcosm experiment, an 
enhancement of both viral abundance and FIC (fre- 
quency of infected cells) in P-enriched samples as a 
result of nutrient stimulation of bacterial growth, which 
in turn enhanced viral activity. However, it is note- 
worthy here that although phosphorus concentration 
was 2-fold higher in Lake Bourget than in Lake Annecy 
(Table 1), no significant difference was recorded in bac- 
terial production between the two lakes (t test, P > 
0.005). Some studies have suggested that nutrient avail- 
ability may have an important influence on viral life 
strategies (e.g. [35,36]). As lysogenic infection is consid- 
ered the most favourable method of bacterial infection 
in water characterized by low bacterial abundance and 
primary production, this may also explain the relatively 
weak stimulation of viral production observed in Lake 
Annecy compared to Lake Bourget [32]. 

In Lake Annecy, and in contrast to viral production, 
the effects of flagellate presence on viral abundance 
seemed to be highly variable between the two periods 
(LAI vs. LA2). This variation revealed viral abundance 
stimulation in early-spring (LAI) and repression in sum- 
mer (LA2), for both treatments (VFA and VF). This 
result could suggest a direct grazing of flagellates on 
viruses during summer. Virivory by flagellates has been 
previously reported [37,38] and according to Domaizon 
et al. [39], all flagellates do not act similarly because of 
large differences between taxon-specifc ingestion rates. 
During our study, heterotrophic flagellates were mainly 
represented by Oikomonas (45 and 48% during LAI and 
LA2, respectively). Also, the grazing impact of flagellates 
on viruses has always been reported to be relatively low, 
resulting in < 4% loss [37,38]. Hence, direct grazing of 
flagellates on viruses was unlikely to explain the repres- 
sion of viral abundance in LA2. Other factors should be 
invoked [36] and would need further investigation. 

Effect of both flagellates and viruses on bacterial activity 

Higher bacterial production in both VF and VFA treat- 
ments than V suggested that grazers and viruses acted 
additively to sustain (directly or indirectly) bacterial 
activity in Lake Annecy and Lake Bourget. Such bacter- 
ial growth and activity stimulation in the presence of 
grazers and viruses has been previously reported in the 
oligotrophic Sep reservoir, Massif Central, France 
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[23,26]. However, the control and reduction of bacterial 
production by the two mortality agents have been 
observed in other aquatic systems [18,21,22]. Such varia- 
bility in possible responses could be due to the initial 
bacterial community composition and environmental 
conditions. 

The increase in bacterial production with the presence 
of both predators (flagellates and viruses) could be 
explained by the fact that grazing activity and viral lysis 
are likely to release inorganic and organic nutrients 
which may stimulate bacterial activity. Obviously, the 
absence of direct measurements of grazing rates of fla- 
gellates on heterotrophic bacteria communities, for 
instance using fluorescently labelled bacteria (FLB) [40], 
prevented us from drawing firm conclusions about the 
grazing pressure of HNF on bacteria and our results 
should be considered in light of that. However, it has 
been suggested that a minimal proportion of 1,000 het- 
erotrophic bacteria for one heterotrophic flagellate is 
characteristic of microbial food webs in which flagellates 
preferentially consume bacteria [39,41,42]. The value for 
this ratio was higher than 1,000 in each treatment (VFA 
vs. VF) and for each experiment (early spring vs. sum- 
mer). Indeed it varied between 1,632 and 3,866 bacteria 
per flagellate in Lake Annecy (mean value: 2,795), and 
between 2,619 and 8,857 in Lake Bourget (mean value: 
5899), suggesting that heterotrophic bacteria were abun- 
dant enough to support the development of the hetero- 
trophic flagellates that were present. 

Seasonal variability in the stimulation of bacterial pro- 
duction seemed to be more important than the trophic 
status variability, with highest mean values recorded in 
summer (+33.5% and +37.5% in Lakes Bourget and 
Annecy, respectively), a period which corresponds to 
low total phosphorus conditions and high temperature 
in surface waters (Table 1). Thus, the input of nutrient 
resources by viral and grazing activities, under such 
summer conditions, is likely to stimulate the bacterial 
community much more than under the cold early-spring 
conditions (temperature = 6-7°C). Moreover, Thomas et 
al. [32] observed that the abundance of HDNA (high 
nucleic acid containing bacteria) is lower in spring than 
in summer in Lake Bourget (less than 40% of the total 
bacterial abundance), and this group is considered to be 
more active in comparison to LDNA (low nucleic acid 
bacteria) [43,44]. This could also explain the low stimu- 
lation of bacterial production in early spring compared 
to that in summer. 

For most experiments (LAI, LB1 and LB2), the stimu- 
lation of bacterial production, at the end of experiments, 
was much higher in VFA than in the VF treatment 
(Figure 4) which could be attributed to an increase in 
substrate availability and regenerated nutrients, resulting 
from grazing pressure of flagellates on both heterotrophic 



bacteria and autotrophic communities, in treatment VFA 
[45,46]. Nevertheless, the opposite situation was observed 
in Lake Annecy during the summer period (LA2 experi- 
ment) suggesting weaker competition by the bacterial 
community for nutrient resources in the presence of 
autotrophs, at this period. As a possible explanation, the 
abundance of autotrophs (represented mainly by picocya- 
nobacteria and PNF) was indeed 2- to 4-fold higher in 
summer than in early spring while bacterial abundance 
was 2-fold lower (Table 1). 

Impact of HNF on bacterial community structure 

We are aware that the DGGE fingerprinting method 
presents some bias and only reflects the microorganism 
populations that are present at relatively high concentra- 
tions. For example, while Muyzer et al. [47] claimed that 
the reported sensitivity of DGGE is 1% of the template 
DNA, Casamayor et al. [48] reported that the number 
of bands is related to the number of populations that 
account for more than 0.3-0.4% of the total cell counts. 
In addition, some other bias such as insufficient or pre- 
ferential disruption of cells during the DNA extraction 
step, amplification bias (chimera and heteroduplex for- 
mation) and band co-migration in the DGGE gel can 
occur and consequently over- or underestimate the 
number of bands. However, such limitations are not 
specific to DGGE and may also be found in other mole- 
cular fingerprinting techniques [49]. Therefore, it must 
be kept in mind that only major changes in the bacterial 
community composition could be monitored using 
DGGE. That is exactly what we observed in this study 
as all sequenced bands belonged to Actinobacteria and 
Proteobacteria, known to be the most dominant phyla 
in lakes [50,51]. Thus our results have to be interpreted 
with caution because the structure of some "non-domi- 
nant" phyla, non-detectable with the DGGE technique, 
could have changed according to the treatments per- 
formed in this study. 

We found that some bands were specific to each treat- 
ment suggesting that some bacterial phylotypes were 
able to develop and thwart the predation pressure. Such 
specificity has already been reported in other experi- 
mental studies [18,21,22]. Phylotypes, observed in both 
VFA and VF treatments, were likely to be resistant to 
both grazing and infection [21,22]. Nevertheless, the 
presence of phylotypes only in VF (not in VFA) might 
indicate sensitivity to the autotrophic activity as a result 
of a weak ability to compete for resources. Phylotypes 
only present when viruses were the exclusive mortality 
agents would probably not be able to deal with the com- 
bined pressure of grazing and viral lysis [21] or were 
strongly susceptible to grazing as already suggested by 
Zhang et al. [22]. Finally, the appearance of bands in 
both VF and VFA treatments could be due to 
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phylotypes benefiting from the presence of predators, e. 
g., via the production of DOM or by the removal of 
competitors. 

A few previous experimental studies have shown that 
viruses can influence the presence or absence of specific 
bacterial phylotypes but, typically, result in a reduction 
of the number of detected phylotypes [53,54]. In con- 
trast, we observed during the summer period an 
increase in the apparent richness when viruses were the 
exclusive mortality agents (i.e. the number of detectable 
bands) giving support to the "killing the winner hypoth- 
esis". The stimulation of bacterial diversity in the pre- 
sence of viruses was also reported in other lacustrine 
systems by Weinbauer et al. [21] and other experimental 
studies performed in coastal marine systems observed 
the same trend [18,22]. However, the relative stability of 
the apparent richness during early spring experiments, 
in treatment V, highlighted the seasonal variability of 
virus effects on bacterial diversity. This high variable 
impact of viruses upon bacterial community structure, 
already reported by Hewson and Fuhrman [54], could 
suggest the influence of stochastic processes. 

Since no decrease in the number of bands was 
observed in either treatment VF or VFA, our result 
could not support the hypothesis of Miki and Yama- 
mura [28] according to whom grazing on infected cells 
"Kills the killer of the winner" and thus reduces bacterial 
species richness. In some cases, the combined effect of 
viruses and flagellates on bacterial fingerprint diversity 
was more consistent than the effect of viruses alone, 
suggesting that both predators acted additively to sus- 
tain apparent richness. According to Zhang et al. [22] 
the 'killing the winner' hypothesis is mediated by both 
predators and not just by one type of predator (viruses). 
Thus, all predators (viruses and flagellates) could act 
additively in controlling the winners of the competition 
for resources and caused an increase in detectable phy- 
lotypes. In addition, stimulation of bacterial production 
and related viral lysis also suggested input of nutrients 
and substrates from grazing and lysis activities which 
may decrease the competition pressure within bacterial 
community, thereby increasing the competitiveness of 
the minor phylotypes [23]. 

The effect of both predators on the bacterial diversity 
was not apparent in all experiments, suggesting more 
variability and complexity in the interactions between 
bacterial diversity, viruses and grazers than hitherto 
assumed. Diverse patterns between predators and bac- 
terial diversity were reported in other studies [18,19,55]. 
Such variability could be explained by the change in the 
balance between bacterial production and protistan 
grazing [56] or to chaotic behaviour due to competition 
among predators for the same prey [28]. Overall, pre- 
vious work performed in both Lakes Annecy and 



Bourget, indicated that the strong complexity of the 
combined physico-chemical and biological parameters 
(with a larger effect of abiotic factors) is mainly respon- 
sible for the evolution of the bacterial community struc- 
ture [57]. 

Conclusion 

Many forms of interaction exist between the various 
components of the microbial loop including the viruses. 
Our results highlight that viruses and flagellated grazers 
are likely to act together, synergistically, on bacterial 
activity and community structure. Even at the commu- 
nity level, interactions between bacterioplankton, viruses 
and grazers are thus much more complex than hitherto 
assumed. More than ever, additional studies are needed 
to fully assess the factors responsible for the variability 
in the interactions between grazers, bacteria and viruses, 
especially in freshwater ecosystems, as well as their eco- 
logical significance for the microbial community struc- 
ture/role and whole ecosystem functioning. 

Methods 

Study sites and sampling 

Water samples were collected from the two largest nat- 
ural lakes in France. For the purpose of this study, 40 
litres of water samples were collected near the surface 
(i.e. 2 m) using a water pump and large tubing on 26 
March and 10 July 2007 in Lake Annecy (referred to 
later as LAI and LA2, respectively) and on 02 April and 
17 July 2007 in Lake Bourget (i.e. LB1 and LB2). In this 
way, for each period, samples were separated by only 
one week between the two lakes. 

Physicochemical variables 

Total organic carbon (TOC) and nutrient concentrations 
(NH 4 , N0 3 , P0 4 , total phosphorus) were measured at 
each station and date, according to the standard French 
protocols AFNOR (details available at http://www.dij on. 
inra.fr/thonon/les_plateaux_techniques/le_laboratoire_- 
de_physico_chimie). A conductivity- temperature-depth 
measuring device (CTD SEABIRD SAB 19 Seacat profi- 
ler) and a Chlorophyll fluorescence Fluoroprobe (BBE 
Moaldenke, Germany) were used to obtain vertical pro- 
files of water temperature, conductivity, dissolved oxy- 
gen concentration and chlorophyll a fluorescence. 

Size fractionation approach 

Immediately after sampling, samples were pre-filtered 
through a 60- [im mesh screen, followed by pre-filtration 
through Nucleopore membranes (< 5-(im pore size) 
under low differential pressure (< 50 mm Hg) in order 
to exclude large eukaryotes. We could thus focus our 
attention on the small eukaryotes, autotrophic and het- 
erotrophic prokaryotes and viruses. A third of the pre- 
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filtered sample was then filtered through 1.6-u.m pore 
size to yield a total free-living bacteria and 'grazer-free' 
containing fraction, which was confirmed by detailed 
microscopic examination at the beginning and at the 
end of the experiments. The remaining pre-filtred sam- 
ple was divided into two parts; one of them was kept in 
a black box (simulating darkness) to inhibit the auto- 
trophic activity. Therefore, three combinations of treat- 
ments were performed: the treatment 'Viruses + 
Bacteria + heterotrophic Flagellates (grazers) + Auto- 
trophs' (fraction < 5 |im, referred to as VFA); the treat- 
ment Viruses + Bacteria + Flagellates (grazers)' (fraction 
< 5 urn put into a black box; VF) and finally the treat- 
ment 'Viruses + Bacteria', i.e. without the flagellates and 
the autotrophic community (fraction < 1.6 |im, referred 
as V). 

Samples so transformed were divided into triplicates 
and poured into 2.5 L Nalgene transparent carboys, 
which had been previously cleaned with 1.2 N HC1 and 
rinsed three times with Milli-Q and filtered lake water. 
All the bottles were incubated in the lake at 2 m depth 
for four complete days. Subsamples were taken from 
each triplicate at day 0, 2 and 4 to assess microbial 
abundances and activities, and at day 0 and 4 for the 
analysis of the bacterial community diversity. 

Flow cytometry (FCM) sample analysis 

We used a FACSCalibur flow cytometer (Becton Dickin- 
son, Franklin Lakes, NJ, U.S.A.) equipped with an air- 
cooled laser providing 15 mW at 488 nm with the stan- 
dard filter set-up. Only a few hours after sampling (less 
than 4 h), one millilitre of water was immediately ana- 
lysed without adding any fixative or dye to analyse the 
picocyanobacterial community dynamics and also to 
check for the absence/presence of prokaryotic (e.g. Syne- 
chococcus) and eukaryotic autotrophic organisms in the 
V treatment. Such unfixed samples, kept in darkness in 
refrigerated boxes and at 4°C for a few hours before 
analysis, revealed no significant changes in cell counts 
while this was not true when using either formaldehyde 
or glutaraldehyde (not shown). Fluorescent microbeads 
(Molecular Probes Inc., Eugene, OR, U.S.A.) of l-(im 
diameter were added to each sample as an internal stan- 
dard. Another 1 ml was fixed and used for bacterial and 
viral counts via FCM, according to the protocol 
described in Personnic et al [25]. Briefly, viruses were 
fixed with glutaraldehyde (0.5% final concentration, 
grade I, Merck) for 30 min in the dark, then diluted in 
0.02 nm filtered TE buffer (0.1 mM Tris-HCL and 1 
mM EDTA, pH 8), and incubated with SYBR Green I 
(at a final 5 x 10" 5 dilution of the commercial stock 
solution; Molecular Probes), for 5 min at ambient tem- 
perature, followed by 10 min at 75°C, and then another 
5 min at room temperature, prior to FCM analysis. 



Heterotrophic bacterial counts were performed on sam- 
ples that had also been fixed with glutaraldehyde (0.5% 
final concentration) for 30 minutes, but the samples 
were then diluted in 0.02 |im filtered deep-lake water 
sample, and incubated with SYBR Green I (10~ 4 dilution 
of the commercial stock solution) for 15 min [25] List- 
mode files were analysed using Cytowin [58]. 

Enumeration of flagellates 

50 ml sub-samples were fixed with glutaraldehyde (1% 
final concentration), stained with primuline [59] and 
collected onto black polycarbonate membranes (0.8-|im 
pore size). For flagellates, slides were prepared within 24 
h after sampling and were stored at -25°C in darkness 
to minimise the loss of autofluorescence [60]. Slides 
were observed at a l,250x magnification using an epi- 
fluorescence microscope (Nikon Eclipse TE200) under 
UV light for heterotrophic nanoflagellates and, under 
blue and green light for pigmented nanoflagellates. 

Bacterial production 

The incorporation of 3 H-leucine was determined follow- 
ing the protocol of Kirchman [61]. For each sample, 5 
sterile eppendorfs (2 ml) received 1 ml of sub-sample. 
Two samples were fixed with formaldehyde (1.6% final 
cone.) to serve as controls. Eppendorfs were inoculated 
with known saturating 3 H-Leu (80 nM final concentra- 
tion, specific activity: 73 Ci.mmol 1 ) and incubated in 
the dark for 2 h. Protein synthesis was stopped by the 
addition of formaldehyde (1.6% final concentration). 
Samples were then filtered through a 25-mm diameter, 
0.22-|im pore size membrane (GTTP). The filters were 
then rinsed twice with 5 ml of trichloroacetic acid 
(TCA, 5% final concentration). The filters were placed 
in scintillation vials, allowed to dry and solubilised with 
1 ml of toluene. After adding 3 ml of the scintillation 
cocktail (Hionic Fluor, Perkin Elmer), the radioactivity 
was counted with a Packard Tricarb Liquid Scintillation 
Analyser 1500. Bacterial production, calculated in 
pmoles l" 1 h" 1 of 3 H-Leucine incorporated into protein, 
was converted in |igC l" 1 h" 1 following Simon and Azam 
[62]: BP (ugC l" 1 h" 1 ) = Leu (mmols Leu L" 1 h" 1 ) x 131.2 
x (%Leu) _1 x (QProtein) x ID); with Qprotein = 0.86 
(ratio of cellular carbon to protein); %Leu = 0.073 (frac- 
tion of leucine in protein). ID = 1 (Isotopic Dilution); 
131.2 = Molecular weight of the leucine. 

Estimation of viral production 

We used the dilution technique of Wilhelm et al [63] in 
order to estimate the viral production throughout the 
experiment at day 0, 2 and 4. 50 ml of sub-samples 
were diluted and mixed with 100 ml of virus-free (0.02- 
|im pore size pre-filtered at day 0 and kept at 4°C) lake 
water, and incubated in dark conditions. Triplicates 
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were made and incubated at in situ temperature in the 
dark. One-ml sub-samples were collected at 0, 3, 6, 12, 
18 and 24 h. Viral production rates were determined 
from first-order regressions of viral abundance versus 
time after correcting for the dilution of the bacterial 
hosts between the samples and the natural community, 
a necessary step to account for the loss of potentially 
infected cells during the filtration. Viral production (VP, 
virus ml" 1 h" 1 ) was calculated as proposed by Hewson 
and Fuhrman [64]: VP = m x (b/B) where m is the 
slope of the regression line, b the concentration of bac- 
teria after dilution, and B the concentration of bacteria 
prior to dilution. 

We estimated the number of lysed bacteria (cell ml" 1 
h" 1 ) during the viral lysis activity by considering an aver- 
age burst size (27) previously estimated for Lake Bourget 
[7,65] with the number of lysed bacteria = Viral produc- 
tion/Burst Size [66]. 

In order to show the effect of the presence of flagel- 
lates on the dynamics and activities of both hetero- 
trophic bacteria and viruses, we calculated the 
stimulation of the different parameters presented above 
(both abundance and production) in treatments VF and 
VFA (as proposed by Bonilla-Findji et al. [18] and 
Zhang et al. [22]). The stimulation corresponds to the 
difference in variation between treatments with flagel- 
lates (VFA or VF treatments) and the treatment without 
flagellates (V treatment) between 0 and 48 h, and 
between 48 h and 96 h, respectively. As an example, the 
equation used to calculate the stimulation of bacterial 
abundance in the VF treatment, between 0 and 48 h, is: 

Stimulation^*) = x W0 ) _ ( x W0 ) (1) 

V BA (V)48 / V BA (V)0 / 

Where BA( V fa)o an d BA( V fa)48 are the abundance of 
bacteria in VFA, at the beginning and after 48 h of incu- 
bation and BA( V )o and BA (V )48 are the abundance of 
bacteria in V, at the beginning and after 48 h of 
incubation. 

Net growth and loss rates of bacteria 

Bacterial net growth rates with bacterial predators (rb, 
d" 1 ) and without predators (r, d" 1 ) were calculated 
from the difference in abundances from day 0 to day 2 
(t = 48 h) and from day 0 and day 4 (t = 96 h), assum- 
ing exponential growth. We used the equations: rb = 
(In Nbt - In NbO)/t and r = (In Nt - In N0)/t; where 
NO and Nt are the bacterial abundances (NbO, Nbt = 
with predators (VFA, VF), NO, Nt = without predators 
(V)) at the beginning and after 48 h or 96 h of incuba- 
tion. The loss rate of bacteria due to grazing activities 
were calculated as the differences between the treat- 
ment with (VFA, VF) and without (V) predators: g = r 
- rb [67]. 



Nucleic acid extraction, PCR and DGGE 

Analysis of the bacterial community structure was 
assessed using Denaturing Gradient Gel Electrophoresis 
(DGGE). Bacteria were harvested from approximately 
250 ml water onto 47-mm diameter, 0.2-|im pore size, 
polycarbonate white membrane filters (Nuclepore) after 
a pre-filtration step through 2-[im pore size polycarbo- 
nate membrane filters (Nuclepore) to eliminate large 
eukaryotes and filamentous cyanobacteria. The filters 
were then stored at -80°C prior to nucleic acid extrac- 
tion. Nucleic acid extraction was performed as described 
in Dorigo et al. [68]. Molecular weight distribution and 
purity of the DNA were assessed by 1% agarose gel elec- 
trophoresis and quantified by both visual comparison 
with molecular weight markers in ethidium bromide 
stained agarose gels (rough estimate) and by optical 
density measurements using NanoDrop ND-1000 Spec- 
trophotometer (Thermo Scientific). Such material was 
then stored at -20°C until PCR amplification. 

PCR reactions were carried out using the Eubacteria- 
specific primer 358-GC [47] and the universal primer 
907 rM [69] which amplify the variable V3 region of the 
16S rRNA gene and yield a DNA fragment of ca. 550 
bp. All PCR amplifications were carried out using about 
30 ng of extracted DNA in a 50 \i\ reaction mix contain- 
ing 10 x Taq reaction buffer (Eurobio), 1.5 mM MgCl 2 , 
120 (iM of each deoxynucleotide, 1 |iM of each primer, 
bovine serum albumin (Sigma, 0.5 mg ml" 1 final concen- 
tration), and 1.25 U Taq DNA polymerase (Eurobluetaq, 
Eurobio). PCR amplification consisted of an initial dena- 
turation step of 94° C for 5 min, followed by 30 cycles of 
denaturation at 94°C for 1 min, annealing at 52°C for 1 
min, and extension at 72°C for 1 min, and a final elon- 
gation step at 72°C for 5 min using a PTC100 thermocy- 
cler (MJ Research). Correct size (ca. 500 bp length) of 
PCR products were determined by 2% agarose gel elec- 
trophoresis with a DNA size standard (Low DNA Mass 
Ladder, GIBCO BRL). 

DGGE analysis was performed on PCR fragments, as 
described in Berdjeb et al. [57] using Ingenyphor U-2® 
(Ingeny international) and by using a 40-80% gradient. 
Since all of the replicates (more than 70) could not be 
placed in the same gel, aliquots of DNA extracts from 
the three replicates of each treatment were pooled, but 
only after we had checked similarity in DGGE patterns 
between replicates for all sampling time points. Digital 
images of the gels were obtained using a Kodak DC290 
camera, and were then saved for further analysis using 
the Microsoft Photo Editor Software. The DGGE band- 
ing patterns were analyzed using the GelCompare II 
software package (Applied Maths, Kortrijk, Belgium) 
and after digitalization of the DGGE gels. Briefly, 
banding patterns were first standardized with a refer- 
ence pattern included in all gels. Each band was 
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described by its position (Y, in pixel on the image file) 
and its relative intensity in the profiles (Pi) which 
could be described as the ratio between the surface of 
the peak (r^) and the sum of the surfaces for all the 
peaks within the profile (N). 

Cloning-sequencing 

From the DGGE gels, the bands of interest were excised, 
placed in sterile water and stored at -20°C. Prior to 
cloning, each excised DGGE band was subjected to a 
freeze-thaw cycle and then centrifuged. DGGE frag- 
ments contained in the supernatant were used as tem- 
plate in a second PCR amplification performed as 
described above. The resulting PCR products were 
cloned with an Invitrogen cloning kit (TOPO TA clon- 
ing) according to the manufacturer's instructions. 
Twelve clones were randomly chosen for each band of 
interest. Each clone was verified by PCR using the com- 
mercial primers M13 and finally sequenced (GATC Bio- 
tech). Sequences were then edited, aligned with 
Genedoc [70] and finally checked for chimeras using 
Bellerophon [71] and the Ribosomal Database Project 
(RDP) [72]. Sequences were finally subjected to BLAST 
and the RDP database to determine the level of similar- 
ity with other 16S rRNA gene sequences available in 
Genbanks. 

Statistical Analysis 

Differences between treatments per experiment, per 
time point were tested for significance using parametric 
analysis of variance (ANOVA) including post hoc test 
analysis (Fisher's protected least significant difference 
test). Testing for normality and homogeneity of variance 
was performed, and data transformation was done when 
required (for all data compared per test). Differences 
were considered significant at P value of < 0.05. We 
compared the difference on the stimulation rate of 
abundance and production of both viral and bacterial 
communities according to the seasons (n = 12) and 
trophic status (n = 24) by using paired t test. 
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